Astronomy & Astrophysics manuscript no. ms5513rev2 


© ESO 2008 


February 5, 2008 





Investigating star formation in the young open cluster NGC 6383" 



E. Paunzen, M. Netopil, K. Zwintz 



O 

o 

(N 

> 
O 



> 
o 



O 



Institut fur Astronomie der Universitat Wien, Tiirkenschanzstr. 17, A- 1180 Wien, Austria 
Received 2006; accepted 2006 

ABSTRACT 

Context. By studying young open clusters, the mechanisms important for star formation over several Myr can be 
examined. For example, accretion rate as a function of rotational velocity can be investigated. Similarly, sequential star 
formation triggered by massive stars with high mass-loss rates can be studied in detail. 

Aims. We identified and characterized probable members of NGC 6383, as well as determined cluster parameters. 
Methods. New Stromgren uvby CCD photometry, obtained by us, is presented. This new data, together with Johnson 
UBV and 2MASS data in the NIR, was used to investigate characteristics of pre- as well as zero age main sequence 
cluster members. 

Results. We present Stromgren uvby CCD photometry for 272 stars in the field of NGC 6383 and derive its reddening, 
E(b — y) = 0.21(4) mag, as well as distance, d = 1-7(3) kpc from the Sun. Several stars with NIR excess and objects 
in the domain of the classical Herbig Ae/Be and T Tauri stars were detected. Two previously known variables were 
identified as rapidly-rotating PMS stars. The field population is clearly separated from the probable members in the 
color-magnitude diagram. 

Conclusions. NGC 6383 is a young open cluster, with an age of less than 4 Myr, undergoing continuous star formation. 
True pre-main sequence members might be found down to absolute magnitudes of +6 mag, with a variety of rotational 
velocities and stellar activities. 
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1. Introduction 

Young open clusters offer the opportunity to investigate 
star formation simultaneously for a significant number of 
stars from high to very low mass. Processes like accretion, 
mass-loss and pulsation in the presence of both local and 
global magnetic fields can be studied (Monin et al. 2006). 
The pre-main sequence (PMS) phase of stellar evolution 
has been extensively investigated in the last decade (James 
et al. 2006). With the detection of magnetic fields (Wade 
et al. 2005) and the application of asteroseismic tools (e.g. 
Ripepi et al. 2006, ; Zwintz, Guenther & Weiss 2006) the 
early phases of the stellar evolution begin to reveal their 
mysteries. 

We investigate the southern young open clus- 
ter NGC 6383 located at cv(2000.0) = 17 h 34 m 48 s , 
(5(2000.0) = -32°34'00"; Z = 355°690, 6 = 0°041 which 
has been a target of several investigations (see Zwintz 
et al. 2005 for a summary). This aggregate is especially 
interesting because it contains several variable stars and, 
due to a very massive binary system in its core, star 
formation is still ongoing (Rauw et al. 2003). 

We present new Stromgren uvby photometry of several 
fields of NGC 6383, together with a detailed analysis of the 
available 2MASS data (Skrutskie et al. 2006). The color- 
magnitude diagram clearly shows continuous star forma- 
tion over about the last 4 Myr. We have derived the cluster 
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parameters and discovered two rapidly-rotating PMS stars. 
These are very interesting objects for further observations. 
An analysis of the NIR data shows that PMS members are 
present to spectral types of M0 or even cooler. Some mem- 
bers show a NIR excess and are located in the photometric 
domain of Herbig Ae/Be and classical T Tauri stars. Based 
on all available information from our observations and the 
literature, a list of bona-fide members of NGC 6383 was 
compiled. 

2. Observations and reduction 

The observations were performed during the night of 
13./14.06.2004 at the 3.6m telescope (ESO-La Silla), 
with the EFOSC2 (fast modus) and the Loral/Lesser 
2048 x 2048 pixel (1 pixel = Of! 157) CCD, which yields a 
5' field-of-view using a standard Stromgren uvby filter set 
(Observer: M. Netopil). 

As NGC 6383 is believed to be rather extended, about 
20', we have chosen five overlapping fields in the cluster area 
to cover the most interesting objects. In total, 33 frames 
in each filter were observed with integration times ranging 
from 5 to 300 seconds. 

The bias-subtraction, dark-correction, flat-fielding and 
point-spread-function fitting were carried out within stan- 
dard IRAF V2.12.2 routines. 

No offsets in instrumental magnitudes between the dif- 
ferent overlapping fields were detected, within the photo- 
metric and transformation error limits. 

The accurate uvby photometry of IC 4651 published 
by Meibom (2000) was used to transform our instrumental 
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Fig. 1. The standard relations of our instrumental magni- 
tudes (suffixes "inst" ) to the highly accurate uvby photom- 
etry of IC 4651 published by Meibom (2000). 
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Fig. 2. The five observed fields in NGC 6383. The sizes 
(by area) of the symbols are inversely proportional to the 
^-magnitudes. 



our internal numbers, the WEBDA and 2MASS identifiers, 
J2000.0 coordinates, the X and Y coordinates within 
our frames, the observed (b — y), mi and C\ values with 
their corresponding errors, V magnitudes, the {B — V) 
colors from the literature and the number of observations, 
respectively. 



magnitudes to standard ones. We have observed IC 4651 at 
different airmasses resulting in 25 individual frames. After 
correcting all frames for the airmass, the following standard 
relations were derived (Fig.[JJ). The errors in the final digits 
of the corresponding quantity are given in parentheses. 



u st = -1.507(3) + 0.985(2)u inst (1) 

v st = -0.441(3) + 1.000(2)tw (2) 

b st = -3.220(3) + 0.986(2)6 <nst (3) 

Vst = -2.262(3) + 0.991 (2)j/ inat (4) 

V st = -2.239(4) + 0.990(2)y mst (5) 



Those relations were applied to the instrumental magni- 
tudes of the stars in the five fields of NGC 6383. 

For a final test, we compared our uvby photometry with 
data from Eggen (1978) for six stars in both our sample and 
his. He used a narrower v filter which only effects mi and 
c\. The values agree very well within the error limits. 

In total, we observed 272 stars, of which 105 could not 
be measured in Stromgren u because of their relative faint- 
ness and the limited quantum efficiency of the instrument. 

The table with all the data for individual cluster 
stars, as well as nonmembers, is available in electronic 
form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5), |http://cd'sweb.u-strasbgfr/Abstract.html| 
or upon request from the first author. The table includes 



3. Results 

Throughout this paper we have adopted the n umbering 
system from WEBDA ( [http://www.univie.ac.at/we bda), 
only in some cases we also list the numbers according to 
FitzGerald et al. (1978, F#) and Zwintz et al. (2005, Z#). 

NGC 6383, which is associated with a large H II region, 
was targeted for several investigations in the past. They 
have been summarized by Zwintz et al. (2005) in some de- 
tail. The basic parameters of this cluster, from different 
sources, are: E(B - V) =0.31(3) mag, 1.4 < d< 2.5 kpc and 
logt<4Myr. Its apparent diameter is still a matter of de- 
bate. FitzGerald et al. (1978) concluded on the basis of star 
counts (see Figure 1 therein) that only a core radius of 2' 
should be considered, whereas Zwintz et al. (2005) give a 
value of 10', which is based on Lynga(1987). Kharchcnko 
et al. (2005) list a core and cluster radius of 4'8, and 15', re- 
spectively. Their age estimate of 5 Myr is compatible with 
the literature whereas their distance of 985 pc seems too 
low. 

Figure^] shows the observed five fields with the location 
of members and non-members. The sizes (by area) of the 
symbols are inversely proportional to the 1^-magnitudes. 
Inspecting the area surrounding NGC 6383 in the Digitized 
Sky Survey (DSS), one can see several empty areas, proba- 
bly caused by dark clouds. Dobashi et al. (2005) have also 
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Fig. 3. Johnson BV photometry by Zwintz et al. (2005) 
in comparison with the mean photometric values from 
WEBDA. The mean A(B — V) = (B — V) w - (B - V) z 
for the intermediate region (0.5 < (B-V) < 1.5) is about 
0.07 mag, whereas for hotter and cooler stars there is a sys- 
tematic trend which can be described by a second-order 
polynomial. 



listed that region in their Atlas and Catalog of Dark Clouds, 
reporting a cloud and clumps extending over 0.27° in this 
area. On the basis of this fact, a density profile using, for 
example, the extensive 2MASS data will evoke erratic re- 
sults for the cluster diameter. However, a cluster radius of 
10' to 15' seems to be within the realm of possibility, since 
the confirmed members are distributed over the whole ob- 
served field of ~ 15' (Fig. [2]). From the comparison of Fig. 
[2] with the DSS images, we conclude that we have covered 
only 60% of the complete cluster area, introducing an addi- 
tional bias. Only further photometric observations can help 
establishing the diameter of NGC 6383. 

All members of NGC 6383 which are later than a spec- 
tral type of AO are still in their PMS phase. Rauw et 
al. (2003) studied the cluster in X-rays using data from 
the XMM-Newton satellite. They detected several X-ray 
sources, which were interpreted as signs of young T Tauri 
stars. FitzGerald et al. (1978) concluded that the central 
bright (V^ = 5.7mag) spectroscopic binary (07V + 07V) 
system, HD 159176 (WEBDA No. 1), triggered star for- 
mation within its surroundings. They derive an age of 
2.8(5) Myr for it. 

No estimation of the cluster metallicity exists in the lit- 
erature. There are neither spectroscopic abundance inves- 
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Fig. 4. Color-magnitude diagram assuming E{b — y) = 
0.21 mag and V — My = 12.0 mag. The zero age main se- 
quence is taken from Philip & Egret (1980) whereas the 
pre-main sequence tracks are from Siess et al. (2000). Non- 
members brighter than +1 mag were selected from the lit- 
erature (e.g. FitzGerald et al. 1978). 



tigations nor a photometric determination available. The 
latter is caused by the fact that only stars hotter than 
AO are on the zero age main sequence and no giant or su- 
pergiant is a member of NGC 6383. Photometric metallic- 
ity calibrations are only available for giants (Hilker 2000) 
or stars later than AO (Karaali et al. 2005). A query in 
WEBDA reveals that there are no other open clusters in 
the galactic vicinity of NGC 6383 with a known metal- 
licity. The galactic distance, Rgc, of NGC 6383 is about 
5.0 to 6.0 kpc depending on its true distance from the Sun 
(i?o = 8.0 kpc). Taking a mean abundance gradient of —0.04 
dexkpc -1 (Cunha & Daflon 2005) for the Milky Way, an 
upper limit for the metallicity of — 0.12dex is evident. Such 
a small difference from the solar value does not alter the 
isochrone significantly (Girardi et al. 2002). However, it has 
to be emphasized that the "intrinsic" range of metallicities 
for a constant Rgc is about ±0.25 dex (Chen et al. 2003), 
but data for open clusters in the direction to the galactic 
center is still very rare. 

3.1. Optical properties 

For the further analysis, we have checked the intrinsic 
consistency of the available Johnson UBV photometry 
included in WEBDA. In total, twelve independent data 
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Table 1. A list of the 44 probable members of NGC 6383 that are still in their PMS phase. These objects were selected 
on the basis of several photometric diagrams. The column "Flag" denotes stars that are in the PMS region in the NIR 
diagram (Fig. [6]). 



No(our) 


No(W EBDA) 


(b - y)o 


M v 


Flag 


No(our) 


No(W EBDA) 


(6 — y)o 


M v 


Flag 


8 


348 


0.823 


4.71 


* 


163 


538 


0.997 


5.43 




13 


222 


1.134 


5.77 




164 


335 


1.027 


5.81 




21 




0.868 


5.18 




169 




0.868 


5.76 




23 




0.872 


5.21 




178 


596 


0.563 


3.62 




26 




0.890 


5.52 


* 


179 


607 


0.520 


3.00 




28 


206 


0.999 


5.81 




180 


611 


0.539 


2.80 




36 


444 


0.461 


2.80 




191 


378 


0.545 


3.54 




44 




1.087 


5.61 


* 


193 


379 


0.738 


4.85 


* 


55 




1.216 


5.84 




197 


382 


0.574 


3.36 


* 


61 


343 


0.894 


5.52 


* 


206 


251 


0.537 


3.39 




73 


199 


1.074 


5.20 


* 


211 


535 


0.728 


4.57 




86 




1.296 


6.07 




217 


388 


0.494 


3.08 


* 


87 




0.828 


4.13 




226 


604 


0.585 


3.09 




90 


543 


0.778 


4.44 




231 


609 


0.748 


4.13 




106 


337 


0.824 


5.41 




273 


128 


0.592 


3.76 




116 


541 


0.639 


3.44 


* 


277 


503 


0.544 


3.06 




118 


546 


0.832 


4.49 


* 


297 


313 


0.494 


2.97 




119 




1.335 


6.00 


* 


301 


571 


0.753 


4.50 




120 




1.175 


6.03 


* 


307 


566 


0.837 


4.17 




126 


134 


0.888 


5.40 


* 


308 


312 


0.497 


3.29 




136 


539 


1.305 


5.75 




312 


310 


0.753 


4.91 




143 




1.328 


6.14 




328 


309 


0.432 


2.56 





sources are available. The only systematic offset was found 
for the (B — V) data by Zwintz et al. (2005). One has to 
keep in mind that the cited photometry was not obtained 
for an astrophysical analysis of cluster properties, but to 
detect pulsating members. Hence, the reduction was tuned 
for high relative photometric accuracy on time scales of 
hours. Not surprisingly, a re-reduction of the data in con- 
formity with Johnson standards arrives at slightly differ- 
ent values. Using the method described in Mermilliod & 
Paunzen (2003), a second-order polynomial fit to the de- 
viations (middle panel of Fig. E]) provides an intrinsically 
more accurate transformation to the Johnson system. We 
have not included these data in our further analysis. 

The PMS evolutionary tracks by Siess et al. (2000) were 
fit to the cluster's color- magnitude diagram, and used to de- 
termine the cluster's distance and reddening. These tracks 
are listed for luminosity and effective temperature, respec- 
tively. The transformation of luminosity into absolute mag- 
nitude is straightforward. More of a problem is the trans- 
formation of effective temperature into Stromgren (b — y) 
because several calibrations are available. First of all, the 
empirical ZAMS listed by Philip & Egret (1980) was used to 
convert the ends of the PMS tracks into the (b — y) and My 
plane. We found that the calibration by Hauck & Kiinzli 
(1996, Table 1) satisfies the color-magnitude diagram best. 
The results by Clem et al. (2004) confirm our choice. The 
final transformation is given as 



G eff = +0.564(8) + 2.961(23)(6-j/) (b - y) < (6) 
Q ef f = +0.564(9) + 0.822(3)(6 - y) (b -y) >0 (7) 



including a small shift of +0.02 mag compared to the orig- 
inal equations to match the standard ZAMS line by Philip 
& Egret (1980). 

Figure 0] shows the color-magnitude diagram with an 
adopted reddening of E(b — y) = 0.21 mag and a dis- 
tance modulus V — My of 12.0 mag. This corresponds to 
E(B — V) = 0.29 mag and a distance of about 1.7 kpc. 
The objects marked with filled circles are members ac- 
cording to the literature, asterisks represent the possible 
PMS members taken from Table [TJ The stars with abso- 
lute magnitudes brighter than +lmag are identified non- 
members taken from the literature (e.g. FitzGerald et al. 
1978). Figure [5] shows the (U - B) versus (B - V) for 
cluster members, with the standard line from Girardi et al. 
(2002), computed for E(B — V) = 0.29 mag and an error 
limit of ±0.05 mag. This error band represents the main 
sequence up to (B — V)q = +0.1 mag very well. Cooler ob- 
jects are still in their PMS phase and can have a peculiar 
behavior in this diagram. The ±0.05 mag error corresponds 
to AE(b~y) =0.036 and AA V = 0.155 mag. In Fig. El two 
outliers are noticeable, the objects No. 17 and 55. These two 
stars show no other abnormal behavior, which prevents any 
reliable conclusion about their nature. The star No. 55 is 
variable according to Zwintz et al. (2005) with an ampli- 
tude of about 25 mmag, which is too small to explain the 
apparent shift. 

The observed color-magnitude diagram (Fig. [5]) shows 
very interesting features. First of all, there is a lack 
of objects with absolute magnitudes between +0.5 and 
+2.0 mag. The only exception is object No. 264, for which 
a large reddening in the NIR is evident (Fig. El)- This 
gap seems to be the borderline between the members of 
NGC 6383 and the field population. The members in the 
PMS phase lie all within the isochrones between 1 to 4 Myr 
down to M v = 2.5 mag (F0). These are PMS A-type ob- 
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Fig. 5. (U — B)o versus {B — V)o diagram for the members 
(hllcd circles) and possible members (asterisks). The stan- 
dard line is taken from Girardi et al. (2002) computed for 
E{B — V) = 0.29 mag and an error limit of ±0.05 mag. The 
only significant outliers are the stars No. 17 and 55 (known 
as variable). 
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Fig. 6. (J — H)q versus (H — Ks)o diagram. The standard 
lines for the main sequence are from Bessell & Brett (1988), 
with the reddening vectors (dotted lines) according to Rieke 
& Lebofsky (1985). The locus of the classical T Tauri stars 
(dashed line) is from Meyer et al. (1997), whereas the region 
of Herbig Ae/Be and emission line F-type objects is taken 
from Hernandez et al. (2005). 



jects mixed together with Herbig Ae/Be stars (Hernandez 
et al. 2005). The classical T Tauri stars have spectral types 
later than GO (Meyer et al. 1997) which corresponds to 
My w 4.3 mag. 

The estimation of distance error is difficult because 
neither a turn-off point nor red giants are present in 
NGC 6383. That is the reason why widely different dis- 
tance estimations (1.4 to 2.5 kpc which corresponds to 
AA/y«l.lmag for a constant reddening) for this aggre- 
gate can be found in the literature. From Fig. [3] we are able 
to make a rough estimation of the distance error by inspect- 
ing the location of members with respect to the ZAMS and 
PMS tracks. Taking into account the error in reddening for 
both parameters we can still fit all features, within obser- 
vational errors, with a shift of ±0.4 mag. The final cluster 
parameters and their errors are: E(B — V) = 0.29(5) mag 
and d= 1.7(3) kpc. 

3.2. Near-Infrared properties 

The 2MASS data base (Skrutskie et al. 2006) with Near- 
Infrared measurements (NIR) was queried for appropriate 
entries in a radius of 10' around the cluster's center, cover- 
ing the five observed fields. The coordinates of these objects 
were then transformed into a rectangular (X,Y) frame and 



compared with the star positions within our images. We 
have only included objects with 2MASS (JHK) measure- 
ments within two arc seconds of the optical source and with 
unambiguous identifications. This resulted in 262 positive 
detections. 

We have used the 2MASS data to search for signs of NIR 
excess and circumstellar dust. The results of the method 
described by Suchkov et al. (2002), analyzing the (V — Kg) 
versus (6 — y)o diagram, are in total agreement with the 
much more accurate conclusions from the (J — H) versus 
(H — Ks)o diagram, shown in Fig. [51 for all bona-fide and 
possible members. 

The standard lines for the main sequence included in 
Fig. [S] are from Bessell & Brett (1988) with the redden- 
ing vectors according to Rieke & Lebofsky (1985). The dif- 
ferences between the 2MASS photometric system and the 
one from Bessell & Brett (1988) are only marginal, in the 
range of a few hundredths of a magnitude (Carpenter 2001, 
Appendix A). The locus of the classical T Tauri stars is 
from Meyer et al. (1997), whereas the region of Herbig 
Ae/Be and emission line F-type objects is taken from 
Hernandez et al. (2005). Notice that the weak-emission T 
Tauri stars can not be separated from main sequence ob- 
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jects (Meyer et al. 1997). We find five objects with a signif- 
icant NIR-excess: 

- Stars with known NIR excess: Nos. 20, 27 (F4) and 264 
(F6) were targets of the investigation by van den Ancker 
et al. (2000). The amount of excess nicely correlates 
with the deviation from the standard relation for normal 
type stars within our diagram (see their Fig. 2). 

— New detections: No. 13 and 32, both show a moderate 
excess, which makes them interesting targets for a spec- 
troscopic investigation. 

The star No. 264 is highly reddened, probably caused 
by circumstellar material. There are also a few probable 
members which show similar but less extreme behavior. 

For the final step, we have investigated all cool stars 
located inside the PMS tracks (Fig. g]) for which 2MASS 
measurements are available. Fourteen of them are placed in 
the domain of Herbig Ae/Be and T Tauri stars, with thirty 
additional objects located in the region of the standard re- 
lation (Fig.[S|). Ten objects have available (U — B) measure- 
ments (Fig. O, placing them well on the standard relation 
given the cluster's reddening. Table [T] lists the 44 objects 
which are good candidate PMS members of NGC 6383, se- 
lected on the basis of the above mentioned criteria. Further 
spectroscopic and kinematic data are required to establish 
their nature and membership. 

3.3. Two rapidly rotating PMS stars 

We would like to comment on two very interesting stars, 
No. 379 (Z64) and 382 (Z71), published by Zwintz et al. 
(2005), for which our new photometry adds important infor- 
mation. Both stars fall well within the PMS domain in the 
NIR (Fig.[S|) and can be regarded as members of NGC 6383. 
This is supported by the reddening-free [mi] values of 0.414 
and 0.525 mag respectively, which are typical for stars later 
than K (Stromgren 1966). The hypothesis that both arc 
highly reddened B-type stars, as formulated by Zwintz et 
al. (2005), can therefore be excluded. The angular velocities 
of both objects are about 16 radians d -1 or approximately 
30 to 40% of their break-up velocity, if we assume that 
the found periods are due to rotation. Such high values are 
not typical for PMS stars (Herbst et al. 2002) but have 
been observed in similar objects (Strassmeier et al. 2005). 
Rapidly-rotating PMS stars are most important for under- 
standing the accretion process combined with the effects of 
stellar magnetic fields (Eisner et al. 2005). 

4. Conclusions 

We have investigated the young open cluster NGC 6383 
with several photometric systems, presenting new 
Stromgren uvby CCD photometry for 272 stars. From this 
data we derive a reddening of E(b — y) = 0.21(4) mag and 
a distance of d = 1.7(3) kpc. An upper age limit of 4Myr 
was determined for NGC 6383. Neither a turn-off point 
nor red giants have been detected. 

In the My versus (b — y)o color-magnitude diagram, 44 
probable PMS members could be traced from A to very 
cool M-type objects. At least 14 of them are in the NIR 
(2MASS) domain of classical T Tauri stars. Appropriate 
isochrones clearly show that members are present from 1 to 
4Myr with a distinct separation from the field population 
at absolute magnitudes between +0.5 and +2.0 mag. 



Five stars with a large NIR excess are unambiguously 
identified in a (J — H)q versus {H — Ks)o diagram. 

We also report the identification of two rapidly-rotating 
PMS stars with angular velocities of approximately 30 to 
40% of their break-up velocity. Such objects are still very 
rare and most important for out understanding of the early 
stages of stellar evolution, in the presence of accretion and 
local magnetic fields. 
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